function [c, ceq] = myconst(x)
    c = 9*x(1,:).^2 + 16*x(2,:).^2 - 144;
    ceq = -1;
end

